Effects of Fe substitution on the electronic, transport, and magnetic properties of 

ZnGa 2 04: A systematic ab-initio study 



O 
O 

(N 

5h , 

< 
(N 



Leonardo Pisani, Tulika Maitra, Roser Valcntf 
Institut fur Theoretische Physik, Universitat Frankfurt, D-6 Frankfurt, Germany 

(Dated: February 5, 2008) 

We present a density functional study of Fe doped into the tetrahedral and octahedral cation 
sites of the wide band gap spinel ZnGa2C>4. We calculate the electronic structure for different 
substitutions and discuss the magnetic and transport properties for each case considering differ- 
ent approximations for the exchange-correlation potential. We show that for certain doped cases, 
significant differences in the predicted behavior are obtained depending on the exchange correla- 
tion potential adopted. Possible applications of the doped systems as magnetic semiconductors are 
outlined. 

PACS numbers: pacs 
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I. INTRODUCTION 

In recent years, spintronics has emerged as one of 
the most studied fields of research in the semiconductor 
physics because of the possibility of exploiting both the 
carrier spin and charge degrees of freedom for the storage 
and transport of information in semiconductor devices^. 
Because of the already existing fabrication technology for 
III-V semiconductors, the doping of these systems with 
magnetic ions (especially the cases of Ga(Mn)As and 
Ga(Mn)N) is being intensively investigated^ and a large 
amount of theoretical and experimental work has been 
done in order to understand the underlying mechanism 
for the ferromagnetism in these doped semiconductors m£ 

Recently, an alternative approach for this phenomenon 
has been considered^ by taking a non-magnetic spinel 
semiconductor as host material and by doping it with 
Fe. The spinel structure (stochiometric formula AB2O4) 
with two types of cation sites (tetrahedral A and octa- 
hedral B) offers new possibilities of obtaining ferromag- 
netic order by doping one or both type of cation sites 
with magnetic ions. First experiments^ on the solid solu- 
tion [ZnGa 2 4 }i-x[Fe 3 A ] x with x=0.05,0.I0 and 0.15 
showed that long-range magnetic order is induced with 
Curie temperatures up to 200K. Similar attempts by 
Krimmel et al& on Feo.76ln2.17S4 revealed, on the con- 
trary, a spin glass state at low temperatures 

Along these lines, we investigate the electronic struc- 
ture of Fe doped ZnGa204 within density functional the- 
ory (DFT) with special emphasis on the magnetic and 
transport properties. Since the description of these prop- 
erties may be affected by the choice of the exchange cor- 
relation (XC) potential considered within DFT, we per- 
form a comparative study of the properties by considering 
two different XC approaches: the Local Density Approx- 
imation (LDA} 7 - and the Generalized Gradient Approxi- 
mation (GGA) 8 . We observe that though the electronic 
structure of the parent compound ZnGa204 remains al- 
most unaffected by the change of the XC potential cho- 
sen, the kind of Fe substitution (e.g. whether in A-site or 
in B-site) shows quite different magnetic and transport 
properties depending on the chosen XC functional. 



We will show that the various doping options allow for 
a variety of interesting behaviors which range from spin- 
polarized metal to the so-called 'transport half-metal'. 
Though both properties are of potential interest for spin- 
tronic purposes, we will discuss that alternative choices 
of substituent or host compound may improve the quality 
of these properties. 

We have organized the paper in the following way. In 
section II we present a study of the volume optimization 
of the parent compound within LDA and GGA, which 
sets the reference frame for the doped cases. In the sub- 
sequent three sections (III, IV and V) we analyse the elec- 
tronic structure properties of Fe-doped ZnGa 2 04 within 
both XC potentials in three different limits, namely, Fc 
doped in the Zn position, Fe doped in the Ga position 
and Fe doped into both Zn and Ga positions respectively. 
Finally in section VI we present the summary of our cal- 
culations with a comparison to the available experimen- 
tal data and conclude with suggestions for some possible 
future experiments. 



II. PARENT COMPOUND 



A. Computational details 



We have performed DFT calculations using the 
full-potential linearized augmented plane-wave code 
WIEN2k 9 . In our calculations, we chose APW+lo as 
the basis set. The atomic sphere radii were chosen to be 
1.8 a.u. for Zn, Ga and Fe and 1.6 a.u. for O. Expansion 
in spherical harmonics for the radial wave functions were 
taken up to I = 10 and the charge densities and potentials 
were represented by spherical harmonics up to L = 6 . 
For Brillouin-zone (BZ) integrations we considered a 60 
k points mesh in the irreducible wedge and the modified 
tetrahedron method was applied^ 
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B. Volume optimization 

ZnGa2C>4 crystallizes in a normal spinel structure 
(space group Fd3m (227)) with the primitive (rhombohe- 
dral) unit cell containing two formula units with Zn, Ga 
and O at the positions (1/8,1/8,1/8), (1/2,1/2,1/2) and 
(u,u,u) respectively (where u=l/4 corresponds to a per- 
fect spinel structure). Zn is sorrounded by a tetrahedral 
environment of oxygens while Ga sits on an octahedral 
position. In order to adopt a reliable approximation for 
the XC potential in the description of the parent com- 
pound, we performed a full volume optimization within 
two widely used parametrizations of the XC functional: 
namely, the LDAi and the GGA£ approach. For each 
selected volume we also considered the relaxation of the 
internal coordinate of oxygen (u parameter) as allowed 
by the symmetry in the space group 227, by using the 
damped Newton dynamics method, until a force value 
smaller than 1 mRyd/a.u. was reached. The energy min- 
imum of the various selected volumes defines then the 
optimal volume. Fig. ^ shows the total electronic energy 
as a function of the volume of the rhombohedral unit cell 
for LDA and GGA (V=4 a 3 where a is the conventional 
cubic cell parameter). 

The optimized lattice constant whose experimental 
valueii is a=8.334 A, is found to be underestimated in 
LDA by about 1% and it is overestimated in GGA by 
about 1.5%. In Table |U we show a comparison of tetra- 
hedral (Zn-O) and octahedral (Ga-O) distances between 
the optimized ones and the experimental ones. Note that 
the experimental values are underestimated by LDA and 
overestimated by GGA. This kind of behaviour is gen- 
erally not unexpected when adopting the LDA or the 
GGA approximations in the description of semiconduc- 
tors. For instance, in a recent ab-initic^ study of the 
spinel ZnAl2 04, the LDA lattice constant underestima- 
tion was found to be of 1.1%. The origin of the LDA 
and GGA shortcoming may be ascribed to the inability 
of both approximations to describe properly the binding 
due to long-range forces like the van der Waals interac- 
tion (apart from the inaccuracy concerning exchange and 
correlation). In contrast, optimization of the u param- 
eter for ZnGa204 converges to 0.3861 within LDA and 
0.3866 within GGA, both in good agreement with the 
experimental value of 0.3867. 

We have also estimated the residual and less important 
discrepancies due to the neglect of the zero point and 
thermal motion of the atoms. Assuming the cubic spinel 
to be an isotropic three dimensional harmonic oscillator 
with a spring strength given by the curvature of the en- 
ergy vs lattice constant function, we estimate the ground 
state length fluctuation to be 0.019 A within LDA. This 
represents a deviation of about 0.2% with respect to the 
classical adiabatic value. The thermal expansion from 
zero to room temperature is expected to produce an ef- 
fect on the lattice parameter less than 0.1%, as inferrable 
from the data of Josties et ali^, therefore these correc- 
tions fall within the inaccuracy of the XC approximations 
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FIG. 1: (Color online) Volume optimization curves for LDA 
and GGA XC-functionals. The energy scale has been made 
arbitrary for plotting reasons and the volume scale measures 
the relative value respect to the experimental one. 



TABLE I: Optimized and experimental tetrahedral and octa- 
hedral distances. 





LDA 


GGA 


exp. 


Zn-O (A) 
Ga-O (A) 


1.943 
1.975 


1.997 
2.018 


1.974 
1.991 



used here. 

From the fit of the energy-volume curves of Fig. 2] 
to the Murnaghan equation of stated we extract a bulk 
modulus for ZnGa 2 4 of 217 GPa within LDA and 146 
GPa within GGA. Structure optimizations for the spinels 
ZnAl204 and ZnGa204 have been performed in the past 
within the framework of the shell modeU^ providing a 
bulk modulus of 273 and 237 GPa respectively. Re- 
cently Levy et al~ have measured the bulk modulus of 
ZnAl204 to be 202 GPa well below the value obtained by 
the shell model (273 GPa). Therefore, assuming that the 
shell model bulk modulus is affected by the same amount 
of discrepancy for both ZnAl204 and ZnGa204 and since 
our calculated values for ZnGa204 are lower than that of 
the shell model, we expect our results to be in a better 
agreement to experiment than the shell model result. 

We conclude from the previous comparison that the 
gradient correction to LDA as implemented in GGA 
doesn't produce a significant improvement upon the LDA 
itself as far as the host compound is concerned. However, 
in the following sections we will see that upon Fe substi- 
tution the LDA and GGA outcomes may turn out to be 
quite different. 



FIG. 2: (Color online) Bandstructure and decomposed den- 
sity of states for ZnGa204 within LDA. The path chosen in 
the Brillouin zone is L=(l/2,l/2,l/2), T = (0,0,0), X = (0,1,0) 
and W= (1/2,1,0) in units of 2tv /a. The DOS is given in units 
of states per eV and per atom in all the figures. 



C. Electronic structure 

In this section we present the density of states and 
bandstructure of the parent compound within the LDA 
and GGA approximations. In Fig. [21 both properties are 
shown only for the LDA approach since no significant 
modification is observed within GGA. These results prove 
to be in good agreement with previous calculations^. 

The spinel ZnGa204 is a direct gap semiconductor with 
a measured band gap of about 4.0 eV—. At the T point 
the band gap is calculated here to be 2.7 eV which re- 
minds us of the renown tendency of LDA to underesti- 
mate the electronic band gap. No significant improve- 
ment on the gap value is detected within GGA. The den- 
sity of states (DOS) shows the valence band to be mainly 
of O character with a peaked Zn weight due to the full d 
orbital shell. The Ga d states are confined below -10 eV 
(not shown) and strongly atomic-like. 

It is important here to point out the dispersive feature 
of the band structure at the T point above the Fermi level 
which has mainly Ga and O character. The associated 
effective mass is calculated to be 0.28 times the electron 
rest mass and is rather isotropic; therefore doping with a 
magnetic ion may open the possibility for a high mobility 
electronic current which simultaneously might be fully 
spin polarized, as we will discuss in the next paragraphs. 
Finally, recalling that the DOS weights are given within 
the atomic sphere radii considered in the APW+lo basis, 
the smallness of the weight of the dispersive s band at T 
(Fig. [21 right panel) is due to the large extension of the 
4s state outside the sphere (interstitial region). 



III. Fe DOPED INTO A POSITION 

Substitution of Fe in one of the two Zn sites of the 
rhombohedral unit cell removes the inversion symmetry 
from the space group Fd3m and reduces it to the maxi- 
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FIG. 3: (Color online) Bandstructure and density of states for 
Fe doped into the A-site of ZnGa204 within LSDA. Upper 
panels: majority states (spin up). Lower panels: minority 
states (spin down). Recall that the DOS weights are given 
within the atomic sphere radii. Therefore the small s weight 
at Ef in the DOS panel is due to the large extension of the 
4s state outside the sphere (interstitial region). 



mal subgroup F43m (space group 216). In addition, the 
original 8 equivalent oxygen positions are now split in two 
non-equivalent sets (Wyckoff position 16e)A&. Defining 
the doping concentration with respect to the total num- 
ber of tetrahedral sites, this new structure corresponds 
to a 50% Fe doping. In a tetrahedral crystal field, the 
Fe d levels are split into energetically lower e g states and 
upper t2 g states. 

We have calculated the electronic structure for this sys- 
tem within the spin polarized versions of the LDA and 
GGA approximations. Both approaches lead to very sim- 
ilar results, therefore we discuss here the LDA results 
and comment on quantitative differences with respect to 
GGA. 

In Fig. |3 the density of states and bandstructure are 
shown for the spin-polarized LDA case (LSDA). Since no 
major modification with respect to the parent compound 
is observed as far as the Zn and Ga weights are concerned 
and the size of the band gap (2.7 eV) remains unchanged, 
we show in the DOS only the Fe and O weights (Fig. [3| 
right panels). 

Note the full inclusion of the iron antibonding 3d ma- 
jority (spin up) states within the band gap of the parent 
compound (Fig. [3] upper right panel) , while the 3c? mi- 
nority (spin down) states hybridize with the highly dis- 
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persive s band at the T point (Fig. flower right panel). 
In a previous work by Nonaka et ali^ the authors did a 
molecular orbital calculation of one Fe atom doped into 
the Zn position of a ZnGa2C>4 cluster and found the Fe 
d-states to be located in the band gap of the host cluster. 

Due to the absence of any d character at the Fermi level 
for the majority states, we may therefore argue that the 
conduction properties of this compound are strongly spin 
dependent. In fact, the s conduction in the spin down 
channel may be hindered by the presence of the relatively 
localized d states which are instead absent in the spin 
up channel. The situation partially resembles the case 
of transition metal ferromagnets like Iron, Cobalt and 
Nickel, where the Fermi level lies in the d| band but the 
effective spin polarization in transport is positive, due to 
the more mobile 4s(|) electrons. 

The exchange splitting is comparable in size with the 
band gap (2.7eV) and much larger than the crystal field 
splitting (0.8 eV). The oxygen hybridization is more pro- 
nounced in the majority spin bands than in the minority 
ones, as indicated by the larger weight of O in the major- 
ity states valence band. This is ascribed to the ti g sym- 
metry of the filled majority states which in a tetrahedral 
ligand field have a cr-bonding with oxygen, in contrast to 
the minority ti g states which are here almost empty. 

According to an ionic picture, iron in A position is ex- 
pected to assume a 3d 6 (Fe 2+ ) configuration in a high 
spin state with a magnetic moment of 4 [Ib- The cal- 
culated total magnetic moment per unit cell is 4.01 [Ib 
and is distributed mainly among Fe, 01 and interstitial 
contributions in a proportion 3.19, 0.08 and 0.47 /is, re- 
spectively. Recalling that the iron muffin tin radius con- 
sidered in the calculation is 0.95 A, we expect the inter- 
stitial magnetic moment to be mainly due to the d orbital 
tails leaking out of the atomic sphere. 

Calculations within GGA give a total magnetic mo- 
ment per unit cell of 4.03 /is, and the Fe, 01 and inter- 
stitial moments are 3.28, 0.07 and 0.46 fiB, respectively. 
Looking at the Fe moment one clearly sees a larger de- 
gree of localization within GGA with respect to LDA. 
To confirm this, we examine the total amount of charge 
in the atomic spheres and we find that it is increased 
in GGA respect to LDA. No other important differences 
have been detected. 

Because of the highly localized nature of iron d states, 
one may argue about the importance of inclusion of or- 
bital on-site correlations on top of a mere LSDA ap- 
proach. We have performed a LDA+U calculation and 
considered two different versions of the double counting 
correction, namely AMF ("around mean field") by Cyzyk 
and SawatzkySS and FLL (full localized limit) by Anisi- 
mov et al. 21 (11=4.5 eV and J=l eV). No relevant modi- 
fication respect to the LSDA case has been noted in both 
cases. In fact, the LDA+U method aims at producing a 
fully orbitally polarized ground state by opening a gap 
between the occupied and non-occupied states. In the 
present case both versions result only in increasing the 
crystal field splitting e g -t2 g of the minority states just 



above the Fermi level pushing the upper 3 t2 g bands fur- 
ther above. 

It is also interesting to note that the ionic radii of Zn 2+ 
and Fe 2+ in a tetrahedral coordination are very similar, 
namely 0.6 and 0.63 A. As a consequence, the global 
(lattice parameter) and local (tetrahedrally coordinating 
oxygen) readjustment of the structure due to the sub- 
stitution is expected not to bring major modification to 
the electronic properties above discussed. This is in fact 
confirmed within both the GGA and LDA calculations 
where the force values acting on the O and Ga atoms 
are around 1 mRyd/a.u. and therefore negligible. This 
also confirms the minor importance of the differences be- 
tween the two approaches when Fe is substituted in a 
tetrahedral site. In the next section we shall see that im- 
portant differences appear when Fe is substituted in an 
octahedral site. 

To determine the ground state magnetic order of this 
systems we performed an antiferromagnetic calculation 
doubling the present unit cell. Comparing the total en- 
ergies of antiferro- and ferromagnetic configurations we 
found that the ferromagnetic state is lower in energy by 
about 61 meV. This state was also predicted to be the 
ground state in this doping limit from a proposed phc- 
nomenological model in an earlier communication by two 
of us2£. 



IV. Fe INTO B POSITION 

Substituting iron into 1 of the 4 Ga sites of the primi- 
tive unit cell of ZnGa204 causes the lattice symmetry to 
lower from cubic (space group 227) to the rhombohedral 
space group R3m (166). In this group, the 8 equivalent 
oxygen positions of the parent group (Wyckoff position 
32e) split into 2c and 6h positions, with the latter being 
the oxygen atoms octahedrally coordinating iron. In an 
octahedral environment, the 3d orbitals of Fe are split 
into energetically lower t2 g and higher e g manifolds. 

For the present case we find a significant difference 
between the LDA and GGA outcomes. In fact we per- 
formed volume optimizations for both approaches and we 
find that LDA produces as a ground state an intermedi- 
ate spin state for Fe, i.e. the unit cell magnetic moment 
is calculated to be 1 /ib- Therefore, within a simplifying 
ionic picture, only the t2 g manifold is populated with 3 
electrons in the up sector and 2 in the down one, resulting 
in a Fe 3+ oxidation state in the low-spin state 1/2. 

On the contrary, GGA describes Fe in a high-spin state 
with a unit cell magnetic moment equal to 5 /is- The 
reason for this strong difference may be connected to the 
element Fe itself. According to LDA, Fe is predicted to be 
non-magnetic and with a fee structure while within GGA 
the correct magnetic bec structure is found 2 ^. Therefore 
we may expect GGA to be more reliable than LDA for 
the B site doping. 

Another reason in favour of GGA is the end compound 
zinc ferrite, ZnFe2 04, where Fe replaces all the Ga atoms. 
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Recent bandstructure calculations^* have shown that the 
zinc ferrite behaves as a metal within LDA and as an an- 
tiferromagnetic insulator with a small gap within GGA. 
The transport gap has been measured to be 0.2 eV— in 
support of the GGA outcome. 

We can understand the strong discrepancy between 
LDA and GGA with the help of the Tableland recalling 
the typical Fe-0 distances for the high-spin state of Fe 
(2.2-2.3 A) and for the low-spin state (1.8-1.9 A). We 
see that the Ga-0 LDA distance is a low-spin distance, 
therefore after substitution and after volume optimiza- 
tion the Fe-0 distance, which increases only by 1-2%, 
will remain of the low-spin state type. On the contrary 
the GGA Ga-0 distance is much closer to the high-spin 
region and thus after volume optimization Fe will be in 
a high-spin state. 

In Fig. 0| we present the GGA bandstructure and 
DOS. As in the previous case, the band gap of the par- 
ent compound remains approximately unchanged (~ 2.7 
eV). The majority spin d states of iron (Fig. 0] upper 
panel) are fully occupied, of which the t2g subband com- 
pletely hybridizes with the oxygen valence band and the 
e ff antibonding states set the Fermi level. Concerning the 
minority states (Fig. flower panel), those antibonding 
with oxygen are completely empty implying that the Fe 
is in an oxidation state 3+. 

From Fig. 0]we find thus a zero-temperature insulator 
with a band gap of 0.4 eV. It is interesting to compare this 
band gap with that of the end compound ZnFe204 where 
within GGA it is calculated to be one order of magnitude 
lessM. This can be ascribed to the direct cation-cation 
bonding in the case of ZnFe204. Since the Fe atoms are 
at the center of edge sharing oxygen octahedra, the t2 ff 
orbitals are able to directly overlap among each other 
producing a widening of the t2 S bands due to hybridiza- 
tion. This causes the minority t2 9 bands just above the 
Fermi level to become closer to with a consequent re- 
duction of the band gap. The spin exchange splitting for 
the present doped system is of about ~2.5 eV, similar to 
the previous case, but the crystal field splitting is around 
~1.8 eV, more than the double as for the previous case. 

The unit cell magnetic moment which is calculated to 
be 5.0 fiB within GGA, is distributed mostly among the 
Fe muffin tin sphere, interstitial region and oxygen muf- 
fin tin sphere in a proportion 3.93 , 0.40 and 0.10 [Ib 
respectively. 

Note that in this case, at the T point the d band in 
the minority sector remains above the Fermi level due to 
the oxidation state of Fe (3+) which leaves the minority 
antibonding d states unpopulated. Moreover, due to the 
larger crystal field effect, only the e g bands hybridize 
with the 4s band at V while the t2 9 states are pushed 
down just above the Fermi level. Therefore, if we were 
to reproduce a situation similar to the A site case (see 
previous section) in which the 4s band should be at the 
Fermi level, we would need a transition metal in a 3+ 
oxidation state with 9 electrons in the d shell, namely 
Zn 3+ , which doesn't exist. 




FIG. 4: (Color online) Bandstructure and density of states 
for Fe in the B-site of ZnGa2 04 within GGA. Upper panels: 
majority (spin up) manifolds. Lower panels: minority (spin 
down) manifolds. The path chosen in the Brillouin zone is V 
= (0,0,0), L = (0,1/2,0), Z= (1/2,1/2,1/2) and F=(l/2,l/2,0) 
in units of the trigonal reciprocal lattice vectors. 



In order to determine the ground state magnetic or- 
der, we compared the total energies of ferromagnetic and 
antiferromagnetic alignments of the iron spin. For this 
purpose we doubled the primitive unit cell into a tetrag- 
onal one containing 4 Zn, 6 Ga, 2 Fe and 16 O atoms. 
We have then considered an antiferromagnetic arrange- 
ment, where the spins of the two non-equivalent Fe atoms 
in the tetragonal unit cell are antiparallel to each other. 
This type of magnetic arrangement gives us planes of fer- 
romagnetically aligned ions with an interplane antiferro- 
magnetic coupling. Comparing the total energy between 
this alignment and the ferromagnetic one, we find the for- 
mer already to be lower in energy by 1 eV with respect 
to the latter. 

In conclusion, we find that doping Fe into Ga site 
in a concentration of 25% results in an antiferromag- 
netic semiconductor at zero temperature which could 
show conducting properties with increasing temperature. 
Namely, the t2 S minority states can be thermally popu- 
lated and produce a spin-polarized current. But, due to 
flatness of the t2 fl bands just above the Fermi level the 
mobility of the carriers would be very low. On the con- 
trary, due to the strong difference in the hybridization 
of e g and t2 g states with oxygen, the e g bandwidth is 
much larger than the t^g one, as clearly seen from the 
bandstructure shown in Fig. 0] Therefore, setting the 
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FIG. 5: (Color online) Crystal structures considered for Fe 
doped in A and B position in a concentration of 25 % (see 
text). In both cases FeA is situated at the corners of the unit 
cell and in case (i) the Fes atoms (Fesi-Fesi/Fesi-Fes2) 
distance is 5.83 A while in case (ii) the Fes atoms distance is 
2.92 A. 

Fermi level into the majority e g manifold would induce 
a half-metallic behaviour with a fully spin-polarized cur- 
rent. This can be achieved by doping for instance with 
Mn instead of Fe. 



V. Fe DOPED INTO A AND B POSITIONS 



The third case that we considered in this study is 
the doping of iron in both sites A and B in a ratio 1:2 
(Fe^Fe^ ) and in a total concentration of 25 %. Accord- 
ingly, we double the rhombohedral primitive unit cell of 
the parent compound into a tetragonal unit cell contain- 
ing 4 Zn, 8 Ga and 16 O. Defining the doping concentra- 
tion with respect to the total number of occupied cation 
sites, 25 % doping of Fe (in both A and B sites) corre- 
sponds to substituting 1 Zn with 1 Fe and 2 Ga with 2 
Fe in the tetragonal unit cell. Two inequivalent struc- 
tures turn out to be consistent with this sort of doping 
(Figs. 03 (i) and (ii)), both belonging to the space group 
Pm (6). In the structure (i) the Fe^ and Fe_B2 atoms 
form zigzag chains with O atoms in between them (e.g. 
FeA-0-Fes2-0-Fe J 4) along the a direction and the Fe^i 
atoms are connected to the Fe^ through O along the c 
direction. The angles FeA-O-Fe^i and Fe^-O-Fe^ are 
equal to 121° and there is no direct path between Fe^i 
and Fes2- In the structure (ii) the two Fes are crystallo- 
graphically equivalent and their distance is much shorter 
(2.92 A) than in the case (i) (5.83 A), therefore direct 
metal-metal bonding is realized in this situation. 

To clarify the role of the on-site correlation typical of 
the transition metal ions for the present case we have per- 
formed calculations within the LDA+U potential, besides 
the GGA one. We will present here only the LDA+U 21 
results since no major differences have been detected 
from the GGA outcome. (The LDA calculation has not 



been taken in consideration because of the failure of LDA 
in describing Fe in the B case.) The values considered for 
the Hubbard and exchange parameters U and J are 4.5 
and 1 eV, respectively. 

In Fig. |Blwe present the spin polarized density of states 
(DOS) and related bandstructure within L(S)DA+U for 
the structure of Fig. [S](i). We have adopted a ferrimag- 
netic alignment of the Fe spins, specifically the Fe^ spins 
are up and Fe# spins are down. Therefore the major- 
ity states of Fe^ are spin up and the majority states of 
Fes are spin down. The upper panel of the Fig. which 
shows the bandstructure and DOS of the spin up species, 
will contain the Fe^ majority states and the Fes minor- 
ity states, while the lower panel of Fig. (spin down 
species) will contain the Fe^ minority states as well as 
the Fes majority states. 

In Fig. [7|we present a simplified schematic energy level 
to show the approximate location of the c?-bands of Fe^ 
and Fes. 

As the two Fe^i and FeB2 bands have rather similar 
features, we show in Fig. (right panel) the sum of the 
two Fes weights. In the energy range from -10 eV to 
about -3 eV the host compound valence band is found 
to have predominant oxygen character. In this range the 
Fe-0 hybridization occurs mainly in the Fes spin down 
bands (see lower panel of Fig. 0) and less importantly in 
the Fe^ states. In the interval from -3 eV to -1 eV the im- 
purity bands are clearly visible: the spin up bands (see 
upper panel) show the Fe^ e g nonbonding and the t2 9 
antibonding states in the order of increasing energy, and 
the spin down bands (lower panel) are mainly Fes e g an- 
tibonding states spanning the same energy range. At the 
Fermi level (see also Fig. [7J) the Fe^ e g spin down band 
(lower panel) , due to the lowering of the point symmetry 
at the A site, splits into an almost fully occupied band 
and an empty e g band at about 1.5 eV above the Fermi 
level. In the upper panel, the Fes ^ig spin up bands are 
found slightly occupied and above them we distinguish 
the empty Fes e g states. In the same energy range we 
find in the lower panel the split e g spin down bands and 
the t2 ff bands of Fe^. The effect of U will be discussed 
in the next paragraph. 

According to an ionic picture, the oxidation states of 
Fe^ and Fes would be 2+ and 3+ respectively. How- 
ever, due to the small overlap of opposite spin bands at 
the Fermi level, the ionic values are slightly changed (in- 
creased for FeA and reduced for Fes). In fact the mag- 
netic moment of the unit cell is found to be -5.97 [Ib-, 
while in the case of an ionic picture it would be -6 hb, 
since the magnetic moment of Fe^ and Fe^'s are respec- 
tively 4 [Ib and -5 /ib (4+2(-5)). A comparison with 
magnetite, which represents the end compound for Fe 
doping in both A and B sites and it has a fully inverted 
spinel structure, i.e. (FeA 3+ )(Fe B 2+ Fe s 3+ )0 4 tells us 
that in the compound (i) the Fe valence states are be- 
ginning to invert, namely the Fe^ increases its oxidation 
state towards 3+ and Fes decreases it to 2.5+. 
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FIG. 6: (Color online) Bandstructure (left panel) and total 
and decomposed density of states (right panel) (color online) 
for the structure of case Fig. (i) within LDA+U with 
U = 4.5 eV and J = 1 eV. Shown are the results for the 
up spin manifold (upper panel of the figure) and the down 
spin manifold (lower panel of the figure). 
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FIG. 7: Schematic picture showing the relative positions 
of FeA and Fes levels in the bandstructure of the doped 
ZnGa2 04 for the crystal structure of Fig. C.B. and V.B. 
denote the conduction and valence bands of the host com- 
pound. 



An important point concerning the bandstructure 
shown in Fig. [5] is the role of the Hubbard U in the 
scmimetal picture. Since the Hubbard U represents the 
on-site Coulomb correlations and no intersite Hubbard 
V has been included in our calculations, no effect on the 
electronic structure at the Fermi level of Fig. H3 is ex- 



pected by varying U since it contains two sets of bands 
originating from two different sites. In fact, we performed 
a further calculation changing the U value of Fe^i and Fes 
and, as expected, only a relative shift respect to the Fermi 
level is found for occupied (downwards shift) and unoc- 
cupied (upwards shift) states leaving the relative position 
of the up and down bands at the Fermi level unchanged. 
The main effect of the Hubbard U is thus to open a gap 
between occupied and unoccupied states and to increase 
the crystal field splitting which is already taken into ac- 
count by an only LDA description. Therefore in the case 
of Fe,B sites, the Hubbard U just adds to the exchange 
interaction constant since the down states are fully oc- 
cupied and the up states are almost completely empty, 
while in the case of Fe^ site due to the presence of Fes 
as second nearest neighbour, which further splits the de- 
generacy of the e g bands, the Hubbard U transforms this 
split into a sizeable gap. In a pure LSDA treatment the 
FeA e s down bands wouldn't show a relevant splitting, 
thus overlapping with Fe^ t2 g up bands and the conduc- 
tion band wouldn't have an almost pure spin polarization 
nature as in the present case. 

Another important effect to take into account is the 
unit cell expansion or contraction due to the replacement 
of some of the Zn and Ga atoms with Fe. In this case we 
cannot describe Fe,4 and Fes as in a 2+ and 3+ oxida- 
tion states, respectively, since the compound is partially 
inverted, and thus we cannot expect minor changes of 
the structure due to the similar ionic radii of Zn 2+ with 
Fe 2+ and Ga 3+ with Fe 3+ , as seen in the previous cases. 
To have an idea of the structural modification, we recall 
that the ionic radius of Fe 3+ in A site is 0.49 A and it is 
0.78 A for Fe 2+ in B site, while the Zn 2+ and Ga 3+ ionic 
radii are respectively 0.6 A and 0.62 A. We shall discuss 
the atomic forces for the present structure at the end of 
this section. 

Finally, we can attempt a rough estimation of the ferri- 
magnetic ordering temperature by exploiting the analogy 
with magnetite Fe304. The latter is known to have the 
highest ferrimagnetic temperature of 858 K. Considering 
the expression k^T a?~4v / 2J abSaSb 25 , where the cou- 
pling constant Jbb between two B sites and between two 
A sites (Jaa) has been neglected with respect to the mag- 
nitude of Jab, one obtains a good estimation of the Neel 
temperature for magnetite when Jab — 19if . In our 
case we assume a simple multiplicative rescaling effect of 
the magnetic ion concentration on the geometric factor 
4V2 and that the exchange constant doesn't depend on 
the absolute spin value of the two Fe. Calculating the 
total energy differences between the ferromagnetic and 
ferrimagnetic alignment of the Fe spins we obtain an ex- 
change constant Jab value of 25 K which by means of the 
above formula provides us with a transition temperature 
of 260 K. 



We now discuss the case Fig. OJii). It is important 
to note that the driving mechanism here is due to the 
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FIG. 8: (Color online) Bandstructure (left panel) and total 
and decomposed density of states (right panel) for the struc- 
ture of case Fig. Q3 (h) within LDA+U with U = 4.5 eV and 
J — 1 eV. Line description as in Fig. |S| 



Fes-Fes bonding while in case (i) Fes are connected to 
Fe^ through oxygen and no direct Fes-Fes bonding is 
realized. 

Fig. [S] shows the total and decomposed densities of 
states within L(S)DA+U of both spin species. The mag- 
netic ground state is found to have antiferromagnetic or- 
der between Fe^ (majority spin f ) and Fes (majority 
spin I). The spin up bands (upper panel) of Fe^ and Fes 
are completely filled and confined in the valence band 
which has a charge gap of about leV with respect to the 
dow spin manyfold. The spin down bands (lower panel) 
are at the Fermi level for both Fe^ and Fes, indicating 
a mixed valence oxidation state. Similar to case (i), the 
point symmetry at the Fe^ site is reduced by the presence 
of next nearest neighbor Fes which leads to the splitting 
of the FeA e g band into two non-degenerate bands. One 
of them is completely empty (the peak just above E^) 
and the other partially filled (the peak at Ep). The re- 
maining charge partially occupies the antibonding spin 
up Fes bands which have t% g symmetry . The system 
thus prefers to assume a partially inverted structure with 
respect to the fully inverse one of magnetite where only 
Fes antibonding states are at the Fermi level. Due to 
the shorter Fes-Fes distance of 2.92 A the metal-metal 
bonding produces a larger bandwidth of Fes states at Es 
as well as a higher electronic population of Fes states in 
(ii) with respect to (i) as is clear by comparison of the 
DOS of the minority Fes states between case (i) and (ii) 



(see Figs. and |S1 upper panel). We may therefore con- 
clude that the degree of inversion in (ii) is larger than 
that in (i) and we may then argue that the driving mech- 
anism in the inversion of magnetite is due to the direct 
cation(B)-cation(B) bonding. 

Moreover, from the bandstructure (Fig. |SJ) it is clear 
that the spin down Fe^ band at E^ (lower panel) is much 
less dispersive than the spin up Fes band (upper panel) 
resulting in a Mott-Hubbard like behavior of the A site. 
Therefore, we find here a particular kind of half-metallic 
behaviour of the system with respect to the conventional 
full spin polarized half-metals. In fact, both spin species 
are here present at the Fermi level but only one of them 
(I t 2g Fes states) is significantly conducting while the 
other (| e g Fe^ states) is strongly localized. This kind 
of system is called transport half-metal as opposed to a 
conventional half metal where the spin polarization (P) is 
100 %. We calculated the value of P for the present case 
according to the definition^ P n = ^ T "l , -^-k , where N 

is the density of states at the Fermi level for the two spin 
species and v the respective Fermi velocity. From the 
calculated band structure and DOS we derive a value of 
Pi equal to 35 % and of P 2 equal to 81 %. Since the 
bulk current is proportional to N(Ep)vp , the P 2 value 
implies that 90.5 % of the current is carried by the spin 
t electrons, (namely the minority spins of Fes) and a 
transport half-metal is clearly realized. 

However, we should recall that the charge carrier con- 
centration in this case is less than that in magnetite and 
therefore the metallicity of this system is poorer than the 
one of magnetite. 

Finally, some considerations on the stability of the pro- 
posed structures should be done at this stage. Due to the 
lack of implementation of atomic forces calculation for an 
LDA+U potential in the WIEN2k code and aware of the 
fact that the inversion issue in spinels is not intrinsically 
due to on-site Hubbard effects, we have estimated the in- 
stability of structures (i) and (ii) respect to inversion by 
evaluating atomic forces within the GGA potential. The 
structures (i) and (ii) considered here are the GGA opti- 
mized ones and have not been relaxed after replacement 
of Zn and Ga with Fe. In Table ITT1 we show the moduli of 
the forces acting on Fe.4, Fes and the oxygen connecting 
them, which are the most significant ones. Note the dif- 
ference in magnitude between cases (i) and (ii) for Fe,4 
and oxygen due to the larger degree of inversion of (ii) 
respect to (i). Comparing the force directions (not shown 
in the Table) we see that the distance between Fe^ and 
oxygen has to become shorter, testifying the tendence of 
both structures to invert. Clearly the structure (ii) is 
significantly more unstable than the structure (i). 
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TABLE II: Comparison of forces on Fe and O atoms in struc- 



(units 


mRy/a.u.) 






Fe A Fe s 





(0 


20.589 6.623 


36.442 


(") 


40.395 9.240 


60.324 



VI. CONCLUSIONS 

In conclusion, we have performed a systematic ab-initio 
study of Fe doped ZnGa2C>4 to investigate the effect of 
Fe substitution on the electronic structure of the parent 
spinel compound ZnGa2C>4 . For an accurate account 
of the magnetic and transport properties we considered 
various exchange correlation functionals and we observed 
that the electronic properties of the doped system de- 
pend significantly on the exchange-correlation functional 
used, especially when Fe is doped on the B cation site 
or on both A and B cation sites of the spinel compound. 
Whereas in the case of Fe doped in the A site both LDA 
and GGA give qualitatively the same results (close to 
an ionic picture), B-site doping leads to different ground 
state spin of Fe within LDA and GGA. During the course 
of this study we also tried to explore the possibilities of 
obtaining a strongly polarized half-metallic groundstate 
which is a desirable property for spintronic applications. 
We observed that Fe doping in A-site could give rise to 
a spin polarized current due to the presence of highly 
dispersive s band in one spin channel. In B-site, Mn in- 
stead of Fe could lead to a half-metallic state by setting 
the Fermi level in a comparatively wider e g band in the 
spin up channel. Doping both A and B sites with Fe 
leads to a so called 'transport half metal' when there is 
a direct Fe-Fe bonding in the B-sites. 

Finally, recent work on [ZnGaaO^i-afFeaG^Ja; for var- 
ious large Fe doping concentrations, showed the exis- 



tence of a ferromagnetic phase up to 200K. Mossbauer 
spectra were measured and provided a quadrupolar split- 
ting of 0.52 mm/s. Assuming a nuclear quadrupole mo- 
ment for 57 Fe of 0.16 barn, as calculated in refi^S within 
the LAPW method, the resulting electric field gradient 
(EFG) is 3.12 V/cm 2 . Alternatively, assuming a nuclear 
quadrupole moment of 0.2 barn, as calculated in ref^S, 
the EFG is 2.5 V/cm 2 . When Fe is substituted into an 
A site no EFG develops due to the spherical local sym- 
metry at the nucleus. Instead, in the case of Fe doped 
into B site the optimized structure provided us with an 
EFG at the Fe site of 2.11 V/cm 2 , which compares fairly 
well with the experimental values. In the doping case of 
Fe into A and B sites we didn't perform a structure opti- 
mization, and, since EFG is quite sensitive to structural 
changes, we obtain quite different values with respect to 
the experimental ones. 

The Mossbauer experiment of Risbud et al. detected 
the presence of only Fe 3+ . This result would discard the 
option of Fe doped in the A site in favor of Fe doped 
on the B site since from our calculations Fe in the A site 
takes up an oxidation state close to 2+ and not to 3+. In 
fact, in ZnGa204 the site preference of Zn 2+ towards the 
tetrahedral sites is much stronger than the Fe one 3 ^ and 
this can explain why Risbud et al. detected the presence 
of only Fe 3+ , which probably occupies predominantly B 
sites. Alternatively, Mg 2+ has a weaker preference for the 
A site than Ga and Fe, therefore replacing Zn 2+ with 
Mg 24 ~ can provide a practically easier way to populate 
A sites with Fe, and MgGa204 would be a better host 
candidate for the properties presented here concerning Fe 
doping on the A site. 
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